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randomness characteristics of renewable energies. As a result, the power generations of distributed gener- 
ators fluctuate frequently. Furthermore, it is difficult to obtain the parameters of microgrids accurately. 
Thus, the study on small signal stability of microgrids is especially important. This paper presents a 
review of exiting small signal stability methods for microgrids. Besides, a new approach based on Sin- 
gular Entropy and Matrix Pencil is proposed. Meanwhile, its validity and effectiveness are confirmed by 
some examples compared with Prony algorithm. Then, the proposed method is also applied into a micro- 
grid with an improved droop control strategy. Finally, a framework based on the proposed method for 
small signal stability, warning and emergency control online in microgrids is available. 
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1. Introduction tainty caused by meteorological factors, it is hard to connect 


renewable energy sources to utility grid directly [7-10]. By inte- 
grating distributed and renewable sources, energy storage devices, 
a variety of loads, data acquisition and supervisory control devices, 
microgrids realize the interface between the distributed renewable 
sources and utility grid [11-13]. 

At present, the researches on microgrids mainly focus on the 
optimal design of system capacity [14], power generation opera- 
tion and control [15-18], system modeling and simulation [19-23], 

s ; : power quality [24,25] and so on, while rarely on the stability 
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: i [26,27]. However, stable and reliable power supply to loads or util- 
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1 Tel.: +86 13567123512. ity grid is the basic function of microgrids. So, study on stability 
2 Tel.: +86 13906508946. of microgrids is very necessary and important. Unlike transmission 


Recently, the shortage of fossil fuels, vulnerability of large grid 
and the power blackouts caused by it have been more and more 
serious [1-3]. Therefore, distributed and renewable energies such 
as wind, solar and so on gained more and more interests [4-6]. 
However, because of their intermittence, randomness and uncer- 
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Nomenclature 
Pi, Qi active power or reactive power output of the 
ith distributed generator (DG), DG;, whose output 
apparent power is S; =P; +jQ;. Especially, Prop (Qrefi) 
is its reference value. Similarly, Pe (Qg) and Sg are the 
ones supplied by utility grid 

P,Q active power or reactive power after digital filter in 

inverter controller 


Pii, Qu active or reactive power demand of the ith load 

kpp, kpi proportion and integration coefficient of DG; for 
active power droop control 

kqp, kgi proportion and integration coefficient of DG; for 
reactive power droop control 

T, XL resistance and inductance of line per kilometer 

wo rated angular frequency of the microgrid 

Eo rated potential of DG 

Vref three-phase reference voltage vector of DG 

Ep potential of DG injected into the microgrid 

Ts sample time of raw data or sampled data for small 
signal stability modes extraction 

y raw data acquisited from microgrid, and y(i) (i= 1, 2, 
..., N) is its ith point, where N is the length of it 

Y formed Hankel matrix by raw data 


Ey, AE,, AE, the kth Singular Entropy, Singular Entropy 
Increment and Differential of Singular Entropy 
Increment, respectively 


n modal order of raw data 

Ài, fi, & the eigenvalue, frequency and damping ratio of ith 
mode 

SNR signal noise ratio 


network or distribution network, microgrids can operate in isolated 
mode or interconnected mode [12]. This feature brings out a serious 
challenge to stability analysis of microgrids. Stability of microgrids 
also includes two parts. One is small signal stability (SSS) in small 
disturbances sense, and the another is the transient stability in large 
disturbances sense. Except for large disturbances, the power out- 
puts of wind turbines, photovoltaic arrays and other renewable 
sources are time-variant, which can be regarded as small distur- 
bances. What is more, they may couple with load fluctuations too. 
Thus, these kinds of small power fluctuations are so frequent that 
they affect the operation status of microgrids constantly. In a word, 
the difficulty and necessity of study on SSS of microgrids is evident. 

This paper reviews the current approaches for SSS of micro- 
grid. Besides, an improved Matrix Pencil combined with Singular 
Entropy is presented and applied to extract the SSS modes of micro- 
grids. The paper is organized as follows. In Section 2, a review of 
exciting and prosiable methods for SSS of microgrids is available. 
In Section 3, the basic mathematical model of Singular Entropy 
Matrix Pencil is proposed. In addition, the microgrid structure and 
an improved control strategy for DGs are provided. As well, some 
simple examples are done to confirm the proposed method. Then, 
the proposed method is applied to the microgrid for its SSS modal 
parameters extraction. Section 4 is the discussion, a framework 
based on the proposed method for the SSS identification, warn- 
ing and emergency control online in microgrids is presented. The 
conclusions are summarized in Section 5. 


2. Study on mall signal stability of microgrids 


Although there are a few approaches for SSS of microgrids so 
far, their disadvantages are obviously. Therefore, a further work 
is strongly demand and a new field needs to be explored. In [26] 


a method for SSS of microgrids in interconnected mode had been 
proposed. A model presented in [27] is suitable for SSS of micro- 
grids in isolated mode. They both featured the eigenvalue analysis 
approach as their theoretical basis and established SSS models in 
interconnected or isolated mode, respectively. However, results of 
eigenvalue approach for a fixed system may deviate far from the 
actual condition, due to the perturbation of system parameters, 
load fluctuations, power scheduling, distribution network distur- 
bances and other causes of small disturbances in the real world. In 
addition, the calculation process of eigenvalues depends on system 
operation status strictly, unfortunately, which is hard to meet. The 
system’s operation status varies frequently and largely, since the 
power outputs of wind turbines and photovoltaic arrays always 
vary in a large range slowly. As a result, it brings tremendous 
challenges to off-line SSS analysis based on eigenvalue analysis 
approach. What is more, conventional matrix eigenvalue approach 
is based on state space matrix obtained by the local linearization 
model of the complex nonlinear system. If we consider sophisti- 
cated control strategies as well as anumber of nonlinear couplings, 
the eigenvalue approach will not be a good solution. Sometimes, the 
linearization model even may not be available. Besides, the inter- 
action among different modes may cause new oscillation modes 
based on normal form theory [28-30]. Modes interaction is abso- 
lutely ignored by eigenvalue approach, because it is just the results 
of first-order Taylor expansion. Especially, its results of SSS arm to 
the whole system. As a result, we cannot know the specific modes 
for each DG or take emergency control strategies corresponding to 
the right DGs, when there are unstable modes in microgrids. 

As mentioned previously, SSS of microgrids may be calculated 
by eigenvalue theory analytically, but there are some serious dis- 
advantages. On the other hand, the SSS modal parameters of 
microgrids can also be identified or extracted by online data 
directly, which can overcome the disadvantage of eigenvalue anal- 
ysis approach. This kind of methods has a wide range of applications 
in transmission network [31], although they also have some draw- 
backs. Common identification methods are Prony algorithm [32], 
Hilbert-Huang Transform (HHT) method [33] and the Matrix Pen- 
cil method [34,35], and so on. Traditional Prony algorithm uses 
complex exponential functions to fit the raw data, so it may pro- 
duce a large number of false modes beside its sensitivity to noise. 
Although some improved iterative Prony algorithms can partially 
inhibit the effect of noise, they also increase the computation 
complexity [36-39]. HHT based modal identification method can 
extract the instantaneous modal information, but it results to a 
high computational complexity. Furthermore, its end effect needs 
to be overcome [40,41]. Matrix Pencil method is a typical signal 
recognition method, too. Recently, it has been applied to power sys- 
tem low frequency oscillation in [42], as it can overcome the noise 
interference and have low computational complexity [34,35]. Tra- 
ditional Matrix Pencil method requires a threshold to obtain the 
modal order, which is depends on whether the ratio between each 
singular value and the largest one bigger than the threshold or not. 
How to set the threshold is affected by the size of raw data and 
the strength of noise, which is hard to be generalized and has no 
appropriate theoretical guidance. 

A lot of relevant analysis showed that the key of the modal 
parameters extraction from raw data for SSS study is how to deter- 
mine modal order accurately. If the modal order is few than its real 
value, it will miss some modes. In contrast, if the modal order is 
more, it will cause over-fit, which results to false modes and reduces 
the accuracy of the modal extraction results. 

In this paper, Singular Entropy is proposed to determine the 
modal order of the raw data. It can not only obtain the modal order 
whether the noise affects the raw data or not, but also can remove 
the false modes and reduce computational complexity automati- 
cally. This method also uses Matrix Pencil method to extract the 
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modal parameters and analyze the SSS of the microgrids based on 
actual operation data in real world. So the SSS information of micro- 
grids can be obtained real-time and online dynamically. It can adapt 
the operation status and small signal modes varying with time in 
microgrids. On the other hand, it can meet the demand for security 
and stability warning and emergency control. 


3. Singular Entropy Matrix Pencil method and its microgrid 
application 


3.1. Mathematical model of Singular Entropy Matrix Pencil 
method 


Firstly, the Singular Entropy is presented to determine the order 
of the signal modes, which can accurately determine the modal 
order regardless of the noise. Then, Matrix Pencil method is used to 
extract the modal parameters of raw data. Combined with Singular 
Entropy and Matrix Pencil, the proposed method can remove the 
false modes and reduce computational complexity, automatically. 
The process can be described as follows. 


3.1.1. Form Hankel matrix structure 


The acquired or simulated raw data y(i) (i=1, 2,...,N) can form 
a Hankel matrix: 
y) y(2) + ¥(L+1) 
y(2) y(3) + YL +2) 


y=] . , . (1) 
y(N) 

where L is the matrix pencil parameter and a proper chosen L 
can suppress noise in y. Take L=N/4 ~ N/3, usually. For Y, we have 
Y=UDV' according to the singular value decomposed (SVD), where 
singular value matrix D is diagonal De R“—)*(*1), Its element cji 
(i=1, 2, ..., m) is the ith singular value of Y, here m=min{N-L, 
L+1} is the minimum dimension of matrix D. Ue R“-))*(N-) is the 
left eigenvector of Y, while Ve R“*1)*(L+1) is the right one. 


yN- L) y(N Zi +1) 


3.1.2. Determine the modal order 

Singular Entropy had been introduced to detect fault in [43]. In 
this paper, we used it to extract the modal order of raw data y. Now, 
define the Singular Entropy Increment: 


D; Oi : 
= a log — i=1,2,...,m (2) 
ja dt 


where the kth order Singular Entropy Ex is defined as: 


k 
E, = So AE, k<m (3) 
i=1 
So the kth order Differential of Singular Entropy Increment is: 
7 AEk = AE, —1 
AE, = k—(k-1) AE, — AEp_1 k=2,3,...,m (4) 


where AE, indicates the Differential of Singular Entropy Incre- 
ment AE;. When the effective signal saturates, the Singular Entropy 
Increment converges to a bounded value rapidly. So, there is an 
inflection point which is corresponding to the signal modal order. 
After that point, the Singular Entropy Increment is small, which 
can be considered as the effect of noise. Due to this advantage of 
Singular Entropy Increment, no matter how much the noise is, the 
increment of Singular Entropy will converge to a bounded value 
with a noticeable jump when the information of effective signal 
tends to saturation. With the help of Differential of Singular Entropy 


Increment, it is easy to extract the order of the corresponded inflec- 
tion point, namely, the modal order n. 


3.1.3. Form anew singular matrix 
After we determined the modal order n, a new matrix D’ can be 
formed by the first nth singular values of D: 


oa 00.. 0 
0 o 0. 0 

Dsj: < Yon, S (5) 
0 0 0- o 
Ow—L-n)xn 


where D’ e R(N-))x", its first n rows are mainly formed by the first 
nth singular values of D and the last N—L—n rows are zeros. The 
formed matrix D’ can eliminate noise effectively, because it has 
removed the singular values caused by noise. 


3.1.4. Obtain modal damping ratio and frequency 
Now, form two matrixes Y4 and Y3: 


Yı =UD'VT (6) 
Y = UD'V} (7) 


where Vj is formed by the 1 to L rows of the right eigenvector 
V’ corresponding to first n singular values; V2 is formed by the 
2 to L+1 rows of the right eigenvector V” corresponding to first 
n singular values. Based on Y; and Y2, a matrix pencil Y> — ÀY; 
can be derived. Modal parameters are included of the eigenval- 
ues of matrix G = Y{Y2, where Yj is the pseudo-inverse matrix 
for Y;. If the number of non-zero eigenvalues for G is M, namely, 
Ai(i=1, 2, ..., M;M<n), the corresponding modal damping ratio 
and frequency of each mode can be obtained by Eqs. (8) and (9): 


a —sgn(Re(Ai)) 
/1 + (Im(In A;/Ts)/Re(In A;/Ts))* 
Re(In A;/Ts) 

j= (9) 


where å; is the ith eigenvalue of ith mode. £; and f; are the damping 
ratio and frequency of the ith mode, respectively. T; is the sample 
time. According to Eq. (8), the system is unstable if any damping 
ratio &; is negative. That is a simple judgement for SSS of microgrid. 


(8) 


3.2. Microgrid structure and control strategy 


The structure diagram of a microgrid for study in this paper is 
shown in Fig. 1. It integrates 10 kW photovoltaic arrays (DG, ), 3 KW 
micro-gas turbine (DG2) and 2 kW fuel cells (DG3). There are some 
local loads at each DG’s terminal, which are Load 1, Load 2 and Load 
3, respectively. DGs are connected to the point of common coupling 
(PCC) by lines, where load 4 attached. Microgrid connected to utility 
grid at PCC simultaneously. The utility grid is a low voltage distribu- 
tion network, whose rated root mean square value of phase-ground 
voltage is 220 V and its rated frequency is 50 Hz. 

Conventional droop control of DG is shown in Fig. 2. According 
to the active and reactive power droop characteristics of DG, its 
decoupling control feature of frequency and potential amplitude 
can be achieved [44,45]. As shown in Fig. 2, the conventional droop 
control strategy can be expressed as: 


w = w0 — kp(P — Pref) 
ar cor! Ave 


where w and E are angular frequency and potential amplitude of 
the inverter, respectively. Meanwhile, wo and Eo are the ones cor- 
respond to the non-load situation. P and Q are the measured output 
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Load1 


Load2 


Linel 1km 


Line2 1km 
C] 


220V 


Load4 Utility Grid 


Load3 


Line3 1km 
E 


Fig. 1. The structure of a microgrid. 


> * 
Pref P Ore Q 


Fig. 2. Traditional droop control strategy of DG. (a) Active droop characteristic 
for frequency control. (b) Reactive droop characteristic for potential amplitude 
control. 


active and reactive power components of the inverter. Pref and Qref 
are the reference ones. kp and ky are the droop coefficients for the 
active and the reactive power. 

The conventional droop control strategy works and coupling 
control can be achieved, when the line is inductive purely. How- 
ever, if we take into account the resistive components, the angular 
frequency of DG will not only relative to the active power, but also 
reactive power [46]. These resistive components may be line resis- 
tance, the equivalent resistance of the filter inductance, parasitic 
resistance of the power electronic components and so on. Sim- 
ilar results can be reached for the potential amplitude [47]. On 
this situation, it is hard to decouple the complex and nonlinear 
system by the conventional droop control strategy as mentioned 
previously. 

From (10), it is can be observed that the traditional droop con- 
trol can be regarded as a Proportional Integral Derivative (PID) type 
controller, or just a proportional controller exactly, whose input is 
the active/reactive power deviations. In order to achieve the decou- 
pling control of the microgrid with resistive components, we add 
an active, reactive power control cross and an integral component 
and the conventional droop controller or proportional controller 
[46,47]. Then, an improved droop control strategy is proposed, as 
shown in Fig. 3, where a, is the angular frequency for the filter 
whose transfer function is G(s)=,/(s +wp). P and Q are the active 


and reactive power after digital filter, respectively. V,¢r is the refer- 
ence voltage output of DG, which is given by 


Vefa =E sin(ot) 
Vep = E sin(wt — 27/3) 
Ve =E sin(wt + 27/3) 


(11) 


In Fig. 3, Ep is the potential of DG injected into the microgrid. 
Prep and Q,eg are reference power outputs of DG;, which is decided 
by the steady operation condition. G;(s) and G(s) are the transfer 
functions of Proportional Integral (PI) regulators. In this paper, we 
take the follow expression: 


Gls)=k „fi 

1(S) = kpp x (12) 
ko 

Go(s) = kqp + = 


According to Fig. 3(a) and Eq. (12), the improved droop control 
strategy can be written as: 


t 
w = Wo + Kkpp(Pref P) H ka | (Pref — P)dt 
0 
t 
+ kl Qa = Q) F ta | (Qref = One 
0 
t 
E = Eo + kqp(Qrer — Q) + kx | (Quer — Q)dt 
0 


t 
= [onten = P) + tx f (Pref = Par 
0 


The important parameters for the microgrid and controller are 
shown in Table 1. Notation that r and x; are the parameters of 
lines, and the length of lines are shown in Fig. 1. The micro- 
grid receives Sg =2 + j0.6 kV A power from utility grid in the given 
steady-state condition according to power flow analysis. All DGs 
take the improved droop control strategy proposed. 


Table 1 

Parameters of the microgrid system. 
Parameters Preft Qheft Pref Qref2 Preg Qref3 Pg Qg 
Values 10kW 0.4kvar 3 kW 0.4 kvar 2 kW 0.3 kvar 2 kW 0.6 kvar 
Parameters Pi Qu Pi2 Qi P13 Qs Pia Qua 
Values 5kW 0.1 kvar 1kW 0.5 kvar 1kW 0.1 kvar 10kW 1 kvar 
Parameters kpp kpi Kap kqi r XL wo Eo 
Values 0.0005 0.0001 0.0001 0 0.243 Q/km 0.83 Q/km 314rad/s 220V 
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Fig. 3. The improved control strategy and structure of DG. (a) Block diagram for control strategy. (b) The control structure for hardware application of DG. 


3.3. Simulation and results 


3.3.1. Some simple examples 

Now, we consider some simple examples, which can confirm 
the validity and effectiveness of the proposed method. If the active 
power P(t) in per unit satisfies the following relationship: 


P(t) = y(t) = yi(t) + ya(t) + y3(t) + e(t) 


where e(t) is white noise. And the sample time is Ts = 0.01 s. Three 
oscillation signals with damping are expressed as: 


(14) 


yilt) = e794" sin(20 x 27t) (15) 
yo(t) = e935" sin(2 x 27t) (16) 
y3(t) = 0.2e~95" sin(0.5 x 27t) (17) 


There are three oscillation signals in P(t), and each signal corre- 
sponds to the two modes in complex frame. So the accurate modal 
order is 2 x 3=6. When we consider the signal with noise, it is sup- 
posed that the white noise sequence meets normal distribution 
€(t)~ N(0, o). Fig. 4 shows the Singular Entropy Increment in case of 
o =0, o =0.025 and o =0.05, where “o” represents the turning point 
corresponding to the modal order of raw data. From Fig. 4, it is 
clear that the Singular Entropy Increment reflects the information 


of modal order accurately (both are 6). When we consider the sig- 
nal without noise, the Singular Entropy Increment is zero after the 
modal order n. When we consider noise, the Singular Entropy Incre- 
ment also has an apparent transition near the modal order n. After 
n, we regard these AE, is caused by noise, which is very small and 
closed to a bounded value, approximately. So we can use Singular 
Entropy Increment to identify the modal order of raw signal. 


3.3.1.1. Case 1: raw data without noise. In this part, we consider the 
case that raw data is not polluted by noise, and give the modal 
parameters extraction results based on Singular Entropy Matrix 
Pencil. Fig. 5(a) gives the Singular Entropy Increment and its dif- 
ferential. 

The modal order result from Singular Entropy Matrix Pencil is 6, 
and the signal noise ratio (SNR) for the fitting is SNR=585.68 dB. 
In the same case, we use the Prony algorithm with 6 orders 
as a comparison. The results show that the signal noise ratio is 
SNR= 151.19 dB. Table 2 and Fig. 6 show the results in detail. 

From Table 2 and Fig. 6, it is easy to find that the singular Entropy 
Matrix Pencil method can identify the modes of original signal well. 
In addition, it is more accurate than Prony algorithm. Especially, 
the phase result of Prony algorithm for mode 1 is far away from the 
accurate value, namely, the bold value as shown in Table 2. 


(a) 


(c) 04 


0.35 


0.3f 
0.254. 


Fig. 4. Relationship between Singular Entropy Increment distribution and the modal order of raw data. (a) Raw data without noise in the time domain. (b) Differentia 


Singular Entropy Increment AE’. (c) Singular Entropy Increment AE. 


ees, 
fab) 
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[e] 
P 


Z. Zeng et al. / Renewable and Sustainable Energy Reviews 15 (2011) 4818-4828 


AE and AE’ 
o 


== 


= AR 
O 6 order 


(b) 0.4 


AE and AE" 
o 


— AE 


S 
O 6 order 


25 30 


4823 


Fig. 5. Singular Entropy Increment and its differential to obtain the modal order. (a) The results for raw data without noise. (b) The results for raw data with noise. 


Table 2 
A comparison of identification results for raw data without noise. 
Frequency (Hz) Damping ratio Amplitude (p.u.) Phase (rad) 
Real Identified Real Identified Real Identified Real Identified 
Prony algorithm 
1 20 20.0001 0.0036 0.00358 1 1.0045 1.5708 0.3142 
2 2 2.00078 0.0278 0.0278 1 1.0035 1.5708 1.4451 
3 0.5 0.5016 0.0793 0.0793 0.2 0.2005 1.5708 1.5394 
Singular Entropy Matrix Pencil 
1 20 20 0.0036 0.00358 1 1 1.5708 1.5708 
2 2 2 0.0278 0.02784 1 1 1.5708 1.5708 
3 0.5 0.5 0.0793 0.07933 0.2 0.2 1.5708 1.5708 


t/s 


Fig. 6. The identification results for raw data without noise. (a) Comparison between raw data and fitting data by Prony algorithm. (b) Error for Prony algorithm (R-F). (c) 
Comparison between raw data and fitting data by Singular Entropy Matrix Pencil method. (d) Error for Singular Entropy Matrix Pencil method (R-F). R, raw data; F, fitting 


data. 
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(b) 


P/p.u. 


(d) : 


P/p.u. 
Oo 


-0.04 


-0.06 


t/s 


Fig. 7. The identification results for raw data with noise. (a) Comparison between raw data and fitting data by Prony algorithm. (b) Error for Prony algorithm (R-F). (c) 
Comparison between raw data and fitting data by Singular Entropy Matrix Pencil method. (d) Error for Singular Entropy Matrix Pencil method (R-F). R, raw data; F, fitting 


data. 


When the raw data has no noise, the extracted modal parame- 
ters can meet some general requirements for simulation, such as 
oscillation modal identification, damping ratio analysis, evaluation 
of the rationality of the microgrid system structure, control equip- 
ment capacity and location planning and configuration, optimal 
design for control strategies and so on. 


3.3.1.2. Case 2: raw data with noise. In this part, we consider the 
situation that the raw data in (14) has noise. We give modal identi- 
fication results based on the Singular Entropy Matrix Pencil method 
and Prony algorithm. It is assumed that the noise meet ¢(t) ~ N(0, 
0.025). Fig. 5(b) shows the changes for Singular Entropy Increment 
and its differential, while the “o” refers to the modal order identi- 
fied. It is easy to see that Singular Entropy Matrix Pencil method 
can get the modal order exactly. 

When o =0.025, we can identify that the modal order of signal 
is 6 by Singular Entropy Matrix Pencil, and the signal noise ratio for 
fitting is SNR=53.51 dB. Prony algorithm with actual modal order 


is used as a comparison, while its SNR is 9.49 dB. A comparison of 
identification results is shown in Table 3 and Fig. 7. Table 3 shows 
that the identified results of Prony algorithm are still poor due to 
its sensitivity to noise even if the correct modal order is used. 

From the comparison of Table 3 and Fig. 7, since Singular 
Entropy Matrix Pencil can exclude false modes, its identifica- 
tion results are very close to the actual signal. From Fig. 5, it 
is clear that the Singular Entropy Increment of the belt in case 
2 become a little worse compared with case 1 because of the 
noise. Although, two methods can extract the modal parameters 
exactly in case 1 without noise. In case 2, the Singular Entropy 
Matrix Pencil method performances well while Prony algorithm 
is failed. 


3.3.2. Applied Singular Entropy Matrix Pencil method into 
microgrid 

Amicrogrid simulation model is established based on the micro- 
grid structure shown in Fig. 1 and the proposed control strategy in 


Table 3 
Acomparison of identification results for raw data with noise. 
Frequency (Hz) Damping rate Amplitude (p.u.) Phase (rad) 
Real Identified Real Identified Real Identified Real Identified 
Prony algorithm 
1 20 20 0.0036 0.0046 1 1.13 1:5741 0.181 
2 2 2 0.0278 0.223 1 1.91 1.571 1.429 
3 0.5 - 0.0793 - 0.2 - 1.571 - 
Singular Entropy Matrix Pencil 
1 20 20 0.0036 0.0036 1 1 1.571 1.568 
2 2 2 0.0278 0.0278 1 0.99 1.571 1.568 
3 0.5 0.5 0.0793 0.0752 0.2 0.19 1.571 1.564 
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Fig. 8. The power outputs of each DG and utility grid. (a) Active power; (b) reactive 
power. 


Fig. 3. The parameters of the microgrid are shown in Table 1. In this 
part, Singular Entropy Matrix Pencil method is applied to analyze 
its SSS. Supposed that the microgrid is in interconnected operated 
mode, and the results in isolated mode can be obtained analogously. 
The small disturbance is set that the DG; has a 5% active power step 
disturbance at 10s. 

Active power signals are chosen for modal extraction because 
they reflect the dynamic of voltage, current and power angle 
directly. The power outputs of each DG are shown in Fig. 8. We 
extract the modal parameters of the system by each DG’s active 
power signal P; (i=1, 2, 3) and utility grid power Pg in the time 
domain 10-18s. The sample time is T;=0.01 s. The identification 
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Table 4 
The main modes of each DG extracted by Singular Entropy Matrix Pencil. 


Active power SNR (dB) Frequency (Hz)/damping ratio 

Pi 168.24 0.318/0.107, 17.681/0.231 

P2 155.00 0.355/0.020, 5.720/0.636, 17.767/0.184 
P3 147.457 0.356/0.016, 5.658/0.621, 16.944/0.209 
Ps 120.242 0.327/0.051, 0.594/0.049, 18.261/0.513 


results are shown in Figs. 9-12. And the main oscillation modes for 
each signal are shown in Table 4. 

According to Table 4, it can be observed that the damping ratio 
of each mode is positive, that is, the system has no positive real 
eigenvalue and the system is stable. It is worth to notice that each 
DG has high frequency oscillation modes whose frequency is larger 
than 5Hz. These modes reflect the rapid transient after system 
disturbance. On the other hand, each DG has low frequency oscilla- 
tion modes whose frequency is less than 5 Hz. Notation that these 
damping ratios corresponding to low frequencies are small. Thus, 
these modes may be more dangerous and easy to be unstable when 
system operation status changes. The DGs have different main oscil- 
lation modes, and that cannot be obtained by eigenvalue analysis. 


4. Discussion 


The stability of microgrids is the basic requirement for a safe and 
reliable distributed power delivery system. Unfortunately, there 
are a few researches on it. Similar to the transmission network, 
the stability for microgrids is divided into small disturbance stabil- 
ity and transient stability too. However, the operation conditions of 
photovoltaic arrays and wind turbines depend on meteorological 
conditions highly. Their stochastic, intermittent and time-variant 
characteristics lead to the power outputs of renewable sources 
change frequently. These can be considered as small disturbances. 
Visibly, small disturbances are so frequent that the SSS is very com- 
plex and important in microgrids. 

Conventional SSS approaches mainly use the local linearized 
state space model of the nonlinear system. Then the SSS results 
can be obtained by eigenvalues analysis. If all the eigenvalues are 
scattered beside left coordinate in the complex plane, the system is 
stable. On the contrary, if there has a positive real part eigenvalue, 
then the system is unstable. However, due to the unique charac- 
teristics of photovoltaic arrays and wind turbines, the linearized 
model is time-variant frequently even hard to get. Therefore, such 
an approach is greatly flawed for SSS of microgrids. In the other 
word, a further work is strongly demand and a new field need to 
be explored for SSS of microgrids. 

In this paper, a new method for SSS of microgrids is pro- 
posed, which can extract the modal parameters from the observed 
raw data. It is verified by some simple examples and a microgrid 
with improved droop control strategy. According to the proposed 
approach, the system is unstable, if there is any mode with negative 
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Fig. 9. Modal parameters extracted using active power of DG1 (P1). (a) Singular Entropy Increment and its differential to get the modal order. (b) Comparison between raw 
data and fitting data by Singular Entropy Matrix Pencil method. (c) Error for Singular Entropy Matrix Pencil method (R-F). R, raw data; F, fitting data. 
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Fig. 10. Modal parameters extracted using active power of DG2 (P2). (a) Singular Entropy Increment and its differential to get the modal order. (b) Comparison between raw 
data and fitting data by Singular Entropy Matrix Pencil method. (c) Error for Singular Entropy Matrix Pencil method (R-F). R, raw data; F, fitting data. 
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Fig. 11. Modal parameters extracted using active power of DG3 (P3). (a) Singular Entropy Increment and its differential to get the modal order. (b) Comparison between raw 


data and fitting data by Singular Entropy Matrix Pencil method. (c) Error for Singular Entropy Matrix Pencil method (R-F). R, raw data; F, fitting data. 
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Fig. 12. Modal parameters extracted using active power of utility grid side (Pg). (a) Singular Entropy Increment and its differential to get the modal order. (b) Comparison 
between raw data and fitting data by Singular Entropy Matrix Pencil method. (c) Error for Singular Entropy Matrix Pencil method (R-F). R, raw data; F, fitting data. 
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Fig. 13. An application framework integrated the proposed approach in microgrids. 


damping which corresponds a positive real part eigenvalue. Thus, 
this approach can be applied to extract the SSS modal information 
real-time and online, warn and emergency control in microgrids. 
The framework is shown in Fig. 13 and the processes can be 
described as follows: 


1. Required active power data be measured in the terminal of DGs 
firstly. Then, the raw data is transmit to microgrid control center 
(MGCC). 

2. After MGCC obtains a certain number of raw data points, modal 
parameters can be extracted by the proposed approach. 

3. According to the modal parameters, if there is any mode with 
negative damping ratio or the damping ratio is less than a pos- 
itive threshold, the microgrid is judged as unstable or near to 
unstable. Thus, some emergency control strategies should be 
taken or warning should be given according to the contingency 
plans. 


5. Conclusions 


In microgrids, the SSS should attract more attention because 
of the unique characteristics of photovoltaic arrays and wind 
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turbines. This paper gives a review on SSS of microgrids and 
proposes a new SSS method based on Singular Entropy and 
Matrix Pencil. Furthermore, a microgrid simulation model with an 
improved control strategy is used to verify it. Some conclusions 
could be summarized as follows: 


Ac 


of 


. SSS of microgrid is a very important issue of microgrids due to 


the unique characteristics of renewable sources, while there are 
rarely researches on it. 


. Singular Entropy Matrix Pencil method uses the raw data 


acquired from microgrids. So it can get its real-time and online 
SSS information dynamically. 


. The proposed method applies Singular Entropy theory to extract 


the actual modal order exactly. It can exclude false modes, 
reduce the computational complexity and improve the accuracy. 
Besides, it combines with Matrix Pencil which can overcome the 
signal noise very well and extract modal parameters accurately. 


. Simulation results have proved the feasibility, validity and accu- 


racy of the presented method. It shows that the proposed 
approach is suitable for SSS of microgrids and the signal noise 
ratio between fitting data and the raw data is high. 


. A framework of the proposed approach for modal parame- 


ters extraction, monitoring, warning, and emergency control in 
microgrids is given. 
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